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Fluctuation-Dissipation Relations (FDRs) are one of 
the corner-stone of statistical mechanics. They offer a 
direct relation between the fluctuations of a system at 
equilibrium and its response to a small external perturba- 
tion. Classical outcome of FDRs are Einstein or Nyquist 
relations, or, more generally, measures of the susceptibil- 
ity, dissipation coefficient or temperature of the system. 
The hypothesis behind the FDRs restrict their applicabil- 
ity to systems that are close to equilibrium. Theoretical 
extrapolation of the FDRs to systems far from equilib- 
rium is currently a very active area of research In 
this context, experimental tests in several glassy systems 
have evidenced violation of FDRs 2]. Furthermore, gen- 
eral identities about fluctuations and dissipation, theo- 
retically derived for time-symmetric out of equilibrium 
systems 0], have been tested in dissipative (non time- 
symmetric) systems like electrical circuit or turbulent 
flow Q ■ Turbulence is actually a very special example of 
far from equilibrium system. Due to its intrinsic dissipa- 
tive nature, an unforced turbulent flow is bound to decay 
to rest. However, in the presence of a permanent forcing, 
a steady state regime can be established, in which forc- 
ing and dissipation equilibrate on average, allowing the 
maintenance of non-zero averaged velocities, with large 
fluctuations covering a wide range of scales. We use mea- 
surements performed in a turbulent von Karman flow to 
show that there is actually a direct link between these 
fluctuations and the mean flow properties, in a way anal- 
ogous to classical FDRs. This approach provides an esti- 
mate of effective statistical temperatures of our turbulent 
flow. 

Theoretical background and definitions.- Describing 
turbulence with tools borrowed from statistical mechan- 
ics is a long-standing dream, starting with Onsager 
Advances in that direction have been recently made for 
flows with symmetries (2D @, axisymmetric Q) using 
tools developed independently by Robert and Somme- 
ria and Miller Q. They consider freely evolving flows 
described by the Euler equation (no forcing and no dis- 
sipation). The Euler equation conserves the energy and, 



for axisymmetric and shear flows, the helicity. In addi- 
tion, owing to the symmetry, there is conservation of a 
local scalar quantity along a velocity line (vorticity in 2D, 
angular momentum for axisymmetry) resulting in a Liou- 
ville theorem and additional global conserved quantities 
as Casimirs of the local scalar quantity. In the Miller- 
Robert-Sommeria theory, the Euler equation develops a 
mixing process leading to a quasi-stationary state on the 
coarse-grained scale. This state is determined by the ini- 
tial conditions and maximizes a mixing entropy under 
conservation of all the inviscid invariants. For forced 
dissipative flows, the strict conservation of the inviscid 
invariants is lost. However, an equilibrium between forc- 
ing and dissipation can establish itself and the system 
can reach a steady state that is a combination of a sta- 
tionary solution of the Euler equation and fluctuations. 
This steady state is selected by forcing and dissipation. 
Some authors [1, @| have proposed to describe this state 
by maximizing a mixing entropy under only particular 
constraints (see for example Hf and I g in Eq. [I} se- 
lected implicitly by forcing and dissipation. This allows 
the derivation of Gibbs states of the system from which 
one derives general identities characterizing the steady 
states, as well as relations between these steady states 
and their fluctuations. 

We apply this approach to the axisymmetric case 
which is relevant to our experimental device. In that 
case, the variables describing the system are the angular 
momentum a, the stream function ip and the rescaled az- 
imuthal vorticity £ = r~ l uje The representation of 

the flow through (a, ip,£) or through the classical velocity 
components (u r ,ue,u z ) in cylindrical coordinates (r,9,z) 
are equivalent since (u r ,0,u z ) = V x (r ipeg) and 
ue — a jr. Furthermore, r~ 1 d r (r~ 1 d r ijj) +r ~ 2 9 2 ?/' = — £. 
In the inviscid, force-free limit, the global quantities con- 
served by the Euler equation are the energy E, the gen- 
eralized helicities Hf and the Casimirs 7 g , given by [7[: 



E = - I ^iprdrdz + — [ —^-rdrdz, 
2 I 2 1 H 
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)rdrdz, I g = g(a)dydz, (1) 



where / and g are arbitrary functions. For forced dissi- 
pative flows, it has been suggested to conserve only the 
energy E and two particular integrals Hf and I g where 
the functions f and g are selected through forcing and 
dissipation 10fl . For Beltrami flows, in which vorticity is 
proportional to velocity everywhere, i.e. u = A V x u, 
the relevant conserved quantities are the energy E and 
the helicity H so that /(er) = er and g = 0]. 

Let us now apply the statistical mechanics approach 
introduced in the previous paragraph to a Beltrami flow. 
The detailed procedure is described in [l, 11|. The mixing 
entropy of the flow S[p] is defined using the probability 
density p(a, £, r) to have a certain couple of values for 
a and £ at each position r. Because of the Beltrami 
hypothesis, the steady state of the flow maximizes S[p] 
at fixed energy E and helicity H . Then, writing SS — 
(55 E — pSH = (where fi^ 1 is a temperature and p an 
"helical potential") and using two different mean field 
approximations, one finds two relations for the averaged 
fields [3 : 



+ H£G = 0, 



(2a) 
(2b) 



as well as two expressions for the probability density 
/9^(cr, r) and /0 CT (£, r), which happen to be Gaussian func- 
tions. The two mean field approximations consist in fix- 
ing independently £ or a to their time average. The 
thermodynamic coefficients p and (3 have been labelled 
accordingly. Considering the first moment of one gets 
an additional relation for the averaged fields: 



(3) 



Finally, considering the second moment of (a, r) and 
/v(£,r), we obtain relations for fluctuations: 



Pa 
Ha r 



(4) 



Comparing Eq. (|2b() with Eq. ([3]) , we see that the four 
thermodynamic coefficients in these equations obey: 



(5) 



t±o_ = fH 
Pa /% ^ 

Dividing Eq. (|2a[) by r and taking its curl, we ob- 
tain a relation between the poloidal components, i.e. the 
components in r-z plane, of the velocity and vorticity: 
= —p^/P^ZJ^. Then, multiplying Eq. (|2b() by r, 
we obtain the following relation for the toroidal, i.e. az- 
imuthal, components: ue= —pa/PaZJe- Finally, using 
Eq. ((H), we get: 

_ Pa 



Pc 



V x u, 



(6) 



and verify that the averaged flow is a Beltrami flow 
with a constant A = —p a /Pa = — Note that this 

flow minimizes the coarse-grained energy at fixed helicity. 
Additionally, combining Eqs. (j2"a|) and (|2b[) . we find that: 

u^4Jvx(w) = A 2 VxM, (7) 

Pa PC 

namely that a spatial variation of the azimuthal vor- 
ticity creates a poloidal velocity with a proportionality 
constant A 2 . This relation can be used to define a sus- 
ceptibility x = 1/ A 2 which happens to be always positive. 

Eqs. Q can easily be recast into: 



U0 = -r- 



= TT 



1 

/V 
x_ 

Pa' 



(8a) 
(8b) 



predicting uniformity of azimuthal velocity and vor- 
ticity fluctuations. Eq. (|8a[) shows that the azimuthal 
velocity fluctuations define an effective statistical tem- 
perature Eq. i|8b|) links the vorticity fluctuations 
to the susceptibility x an d a vortical effective tempera- 
ture 1/fta- This is a formal equivalent of the Einstein 
relation for the Brownian motion. These equations may 
be regarded as formally analogous to FDRs since they 
link fluctuations, susceptibility and temperature. These 
predictions enable the measurements of turbulence effec- 
tive temperatures through fluctuations of ue and ujg in 
a Beltrami flow. Because variances are positive, (3 a and 
/?£ are always positive (since \ > 0), unlike in the 2D 
situation where x can be negative [la ]. In contrast, p a 
and p^ can take positive or negative values, depending 
on the helicity sign. 

The analogy between our predictions and FDRs can 
actually be pushed forward. An other possible way to 
derive Eqs. (j4|) is to introduce, as in classical statisti- 
cal mechanics, the partition function Z^^ describing the 
Beltrami equilibrium state in each mean field approxima- 
tion: 



1 5 2 \ogZ i _ 1 8a 
Pa Sa 
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(9a) 



(9b) 



Formally, the mathematical objects Sa/S^ and 5!; /Scf 
can be seen as response functions. With this point of 
view, Eqs. (J9j) again reflect a formal analogy with FDRs 
since another classical way to write it down is to link the 
fluctuations of a field to its response to a perturbation. 

The challenge is to face such fluctuation relations with 
experimental data. For this, we use a turbulent von 
Karman flow that has already been shown to tend to 
a Beltrami flow at large Reynolds number ■ 

Experimental flow and measurement techniques. - 
Our experimental setup consists of a plcxiglas cylinder 
(radius R c = 100 mm) filled up with water. The fluid is 
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mechanically stirred by a pair of impellers driven by two 
independent motors in exact countcrrotation. The result- 
ing flow belongs to the von Karman class of flows with 
a mean flow divided into two toric cells separated by an 
azimuthal shear layer. We define the Reynolds number 
as Re = UL/v = 2itFR? c /v, where F is the impeller fre- 
quency and v the water viscosity. Rotating the impellers 
from 2 to 8Hz, we can achieve Reynolds numbers from 
125, 000 to 500, 000. Our two models of impellers, TM60 
and TM73, are flat disks of respective radius 0.925 R c and 
0.75 R c , fitted with radial blades of height 0.2 R c and re- 
spective curvature 0.50 R c and 0.925 R c . The inner face 
of the discs are 1.8 R c appart. Different forcings are asso- 
ciated with the convex or concave face of the blades going 
forward, denoted in the sequel by senses (+) and (— ). Ve- 
locity measurements are done with a Stereoscopic Partic- 
ule Image Velocimetry system [ll| provided by DANTEC 
Dynamics. The cylinder is mounted inside a water filled 
square plexiglas container in order to reduce optical de- 
formations. Two digital cameras are aiming at a meridian 
plane of the flow through two perpendicular faces of the 
square container giving a 2D-three components velocity 
field map. Correlation calculations are done on 32 x 32 
pixels 2 windows with 50% overlap. As a result, each ve- 
locity is averaged on a 4.16 x 4.16 mm 2 window over the 
1.5 mm laser sheet thickness. The spatial resolution is 
2.08 mm. A basic measurement is a set of 5000 acqui- 
sitions at a rate of 4 images per second. From this set, 
we compute time-average and fluctuations of the three 
velocity components. 

From now on, all physical quantities and equations will 
be considered in their non-dimensionalized form using 
F and R c . For each of the four forcing configurations, 
we made between five and seven tests that show no Re- 
dependency for these non-dimensionalized quantities. 

Test of the mean flow relations.- Relations (|2ap and 
(|2b[) can be tested by plotting ip and £r 2 with respect to 
cf for each experiment. Two of them are plotted in Fig. [T] 
As in Ref. [lOj, we focus on the flow bulk, i.e. \z\ < 0.5 
and \r\ < 0.5, where the von Karman flow is close to 
a Beltrami flow. Linear dependencies of Eqs. ^ are 
confirmed [lj| and enable estimation of the slopes f} a / ' u a 
and /?£/ (see caption of Fig. [1] for details and Table 
U for averaged results) . Depending on the forcing, the 
two measurements of fJ/fi differ from 1 to 13%, verifying 
Eq. ([5]) and giving a unique mean value of (3/n or A (see 
Table [H line 3). The quality of the test of Eq. © is 
evaluated by — Pa/t^a (line 4). One can see that 

the two 68% confidence intervals of /%/ and P„jn, a 
always overlap. 

Test of the fluctuation relations.- Now, we turn to ex- 
perimental test of fluctuation relations ([8]) and to com- 
plete determination of the four a priori independent co- 
efficients: Pa, an d u a - Because of Eq. jS]), only 
three of them are independent and the previous test al- 
ready provided a measurement of So, we can use 
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FIG. 1: tp (top) and £r (bottom) versus a for two experimen- 
tal von Karman flow with TM60 impellers at F = 6Hz, sense 
(+) at left, ( — ) at right. Black dots correspond to flow bulk 
data (| 2 1 < 0.5, |r| < 0.5) and define mostly linear functions. 
The slopes of the dot-dashed lines is given by the first or- 
der coefficient of an odd cubic fit of the data. Corresponding 
values of j3 a / ' [i a and are used to compute Table [I] 
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FIG. 2: Evolution of angular velocity fluctuations (top) and 
angular momentum fluctuations (bottom) with the radial co- 
ordinate r for the same flows as in Fig. Q] Black dots corre- 
spond to flow bulk data (| z \ < 0.5, \r\ < 0.5). The correspond- 
ing mean values and standard deviations at each z are plotted 
by red lines and errorbars. Horizontal dot-dashed lines show 
the averages over the flow bulk, i.e. measured values for 
and Ax/74- 

the fluctuation relations of Eqs. {8} to compute the re- 
maining parameters 1//% and /3 CT //i 2 . 

Fig. [2] presents the analysis of the fluctuation rela- 
tions for angular velocity ug and vorticity ujg. On the 
top plots, we test relation (fgaj) . Over the whole flow, 
the velocity fluctuations are roughly constant. The rel- 
ative scattering, which mostly tracks the z-dependance, 
increases with r. Focusing again on the flow bulk, we ob- 
serve a reduced scattering and measure 1 with a sim- 
ple average. On the bottom plots, we perform the same 
analysis for vorticity fluctuations and relation (|8b[) . The 
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Impellers 


TM73 


TM60 


Sense 


(+) 


(-) 


(+) 


(-) 




4.64 ±0.25 


-4.92 ±0.12 


3.76 ±0.28 


-4.11 ±0.31 


Pa/^a 


4.31 ±0.20 


-4.88 ±0.17 


3.55 ±0.20 


-3.61 ±0.23 


< P/fJ, > 


4.47 ±0.22 


-4.90 ±0.15 


3.66 ±0.24 


-3.86 ±0.27 


Pt/Vi- Pal [la 


0.33 ±0.45 


-0.04 ±0.29 


0.21 ±0.48 


-0.50 ±0.54 




0.0452 ± 0.0040 


0.0673 ± 0.0035 


0.0481 ± 0.0056 


0.0922 ± 0.0086 


P.I A 


9.1 ±0.5 


13.4 ±0.7 


10.4 ± 1.4 


16.4 ±0.9 


Pi 


22.1 ± 2.0 


14.9 ±0.8 


20.8 ± 2.4 


10.8 ±1.0 


P. 


2.04 ±0.31 


1.77 ±0.21 


1.22 ±0.31 


0.79 ±0.14 




4.77 ± 0.68 


-3.02 ±0.23 


5.53 ±1.06 


-2.64 ±0.45 


/J-a 


0.47 ± 0.05 


-0.36 ±0.03 


0.34 ± 0.07 


-0.22 ±0.03 



TABLE I: Non-dimensionalized statistical 
coefficients measured in our experiments 
following Eqs. (J5J) and (|5J for the four forc- 
ings studied (two impellers, two senses). 
Raw measurements (lines 1-2 and 5-6) al- 
low to calculate a Beltrami factor (line 3- 
4) and statistical coefficients (lines 7-10). 
Errors are calculated with standard 68% 
confidence intervals as the sum of the er- 
ror on data fits (see Figs. [T] and and 
of the statistical dispersion over the differ- 
ent runs. Each measurement is an average 
over 5 to 7 experiments performed at high 
Re between 1.2 and 5 x 10 5 . 



vorticity fluctuations displays a scattering of the same 
order of magnitude than angular velocity fluctuations. 
Once again, an average over the bulk allows to measure 
Pa/fj,^. We conclude that, for these von Karman flow, 
the variances of azimuthal velocity and vorticity are con- 
stant in the flow bulk. This result is compatible with the 
mean field analysis drawn in the first part of this Letter 
for a Beltrami flow. Combining the two sets of results 
described above, we can independently evaluate (3 a ) 
and /%). Standard error estimates are typically 10%, 
ranging from 5 to 20%. Tablc[T]shows that these two cou- 
ples are not equal to each other and even differ by one 
order of magnitude. This result is discussed in the next 
section. 

Conclusion.- We have derived predictions concerning 
the mean flow and the fluctuations in an axisymmetric 
Euler-Bcltrami flow using tools borrowed from statistical 
mechanics. These predictions, and especially the unifor- 
mity of the variances of ug and ujg (Eqs. [8]), have then 
been experimentally confirmed in the specific case of a 
turbulent von Karman flow at large Reynolds number. 
This result is a priori unexpected since the theory, issued 
from the Euler equation, does not explicitly take into 
account the forcing and the dissipation which implicitly 
determine the form of the steady state. This is an addi- 
tional confirmation that out of equilibrium steady states 
of a real turbulent flow may be described as equilibrium 
states of the Euler equation as suggested in [lfj. Addi- 
tionally, these relations provide two different values for 
an effective statistical temperature 1 //3 of our system de- 
pending on the selected variable. Such non-uniqueness of 
statistical temperature in out of equilibrium systems has 
already been encountered in the context of glassy systems 
[T2I ] . Turbulent flows may be another example of this out 
of equilibrium property. Meanwhile, one could possibly 
use correlation with hydrodynamics properties (such as 
variation with forcing, fluctuation rate or Reynolds num- 
ber to decide whether one of the temperatures we 
measure is more relevant than the other to stand for a sta- 
tistical temperature of the turbulence. Finally, in order 



to test further the physical relevance of the FDR anal- 
ogy, one should perform direct measurements of response 
functions as suggested by the r.h.s. of Eqs. ©■ This is 
a yet unresolved and challenging experimental issue. 
This work was supported by ANR TSF NT05-1-41492. 
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